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We present the equation of state (EOS) of quark-diquark matter in the quark mass- 
dcnsity-dcpcndcnt model. The region of the 2 — D parameter space inside which this 
quark-diquark matter is stable against diquark disassembling and hadronization is deter- 
mined. Motivated by observational data suggesting a high compactness of some neutron 
stars (NS) we present models based on the present EOS and compare the results with 
those obtained with previous works addressing a quark-diquark composition based on 
the MIT Bag model. We show that very compact self-bound stars (yet having > 1M@) 
are allowed by our EOS even if the diquark itself is unbound. 

1. Introduction 

The spectacular advances of astronomical instrumentation in the last decade has 
stimulated a great deal of activity on "old" questions of relativistic astrophysics. 
Among the latter we find the internal structure of compact stars, which is not only 
interesting per se, but also important for a through knowledge of the QCD diagram 
in the low-T , high fi region. 

One of the most simple forms of (indirectly) investigating the nature of high- 
density matter, more precisely of the EOS is by obtaining accurate determinations of 
masses and radii. In fact, leaving aside the possibility of differences due to extreme 
rotation, and adopting General Relativity as a framework, a comparison of the static 
models generated by integration of the Tolman-Oppenheimer-Volkoff equations with 
observed data should reveal a great deal of information about the equation of state. 
Only recently we have achieved sufficient accuracy to attempt some key tests on 
selected objects, although it should be acknowledged that three decades of compact 
star astrophysics had produced quite clever arguments to determine masses and 
radii, even in the cases in which they proved to be wrong. 

A paradigmatic example of the above assertions is the celebrated determina- 
tion of the binary pulsar mass PSR 1913+16, believed to be accurate to the fourth 



*Email: glugones@iagusp.usp.br 



Quark- diquark matter equation of state 



decimal place Id. Other methods based on combinations of spectroscopic and pho- 
tometric techniques have been recently perfectioned, and have confirmed that at 
least one X-ray source is significantly above the "canonical!! 1.44 M Q ; namely Vela 
X-l for which a value of 1.87^0 17 M@ has been obtained Id. Some works on com- 
pact stars have attracted the attention recently because not only the masses have 
been determined, but also indications of the radii are available, suggesting a very 
compact structure. 

On the other hand, from a particle physics point of view it is fairly well estab- 
lished that at high temperatures and/or high baryon number densities confinement 
is effectively disabled. The old question of whether the deconfined matter takes 
place inside "neutron stars" is not yet settled, although a great deal of interest in 
this issue emerged from recent calculation of superconducting QCD phases (two- 
flavor superconducting (2SC) and "color flavor locking" (CFL), see □ for a review). 
In fact, there is growing evidence that paired states or diquarks (e. g. a spin-0, 
color-antitriplet bound state of two quarks) might occur as a component in the 
QCD plasma. This is important for understanding hadron structure and high en- 
ergy particle reactions. Of course, at extremely high densities diquarks are expected 
to lose their identity and dissolve into quarks. In the intermediate phase diquarks 
would be expected to be favored. 

Diquark correlations were originally thought to arise in part from spin-dependent 
interactions between two quarks. Later the non-perturbative interactions were 
claimed to produce a gap substantially greater than in the '80s, possibly as large 
as ~ 100 MeV. Although there is no consensus about diquark details, it seems 
certain that a (ud) quark pair experiences some attraction when in total spin-0 and 
colour-3* states. Many works assume that the diquark correlation is similar to the 
one experienced by electrons in a superconductor, with diquarks being the Cooper 
pairs of QCD □. However, we shall not be concerned with the details of the binding 
of the diquark itself in this work. This is because for compact star astrophysics it 
is perhaps more important to look for the full free energy (per baryon) behavior in 
bulk rather than the mass of the diquark individually. 

2. The equation of state 

The quark mass-density-dependent model (QMDD model) is a phenomenological 
description that attempts to incorporate the confinement and asymptotic freedom 
of quarks as a function of the baryonic densitv of matter. According to lattice 
calculations c3 and string model investigations c3, the quark-quark interaction is 
proportional to the distance so that m oc rig 1 ^ 3 may be assumed LL3. Therefore, it is 
reasonable to adopt the following paramctrizations of the particle masses. For the 
quarks u and d we define as usual 



C 

m u =md = —, 

n B 



(1) 
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For diquarks, we introduce here the mass parameter moo and write the diquark 
mass mo as 

2C 

m D = m Da + m u + m d = m Da + -— ^. (2) 

n B 

The non-perturbative character of the interactions between quarks and diquarks 
is included through this particular mass parametrization with the baryon number 
density. 

Before describing the equation of state in detail we would like to make a num- 
ber of qualitative remarks: To make contact with current calculations we must 
determine to which extent the model describes the expected diquark physics. The 
diquark pairing gap that arises from current calculations □ clearly depends on the 
chemical potential (or on the baryon number density) . Qualitatively the gap is zero 
at low densities, reaches a maximum value in the 10-100 MeV range at intermediate 
densities and goes to zero at high densities. To mimic this behaviour the parameter 
todo should have the same shape (but with negative sign) as a function of, for ex- 
ample, the baryon number density, even though in QMDD models it is introduced 
to simulate confinement only. We shall consider only constant values for tudo m 
this work. As we shall see below this approximation is a reasonable starting point 
for the physics of interest to us here. On the other hand we must note that recent 
lattice simulations show that the diquark could be unbound El The latter situation 
can be described by the present model with positive values of m^o, which are also 
included in the analysis (see below). We must remark that for negative values of 
moo the diquark mass becomes negative (as expected) for some baryon number 
density beyond which the model is no longer valid. However, this ng is very high 
for typical values of moo (for example, moo — -100 MeV gives nB,Hm greater than 
12 x no, with uq — 0.16 fm~ 3 the nuclear saturation density. We note that our 
quark-diquark phase may be directly associated to the phase explicitely addressed 
by Rapp et al. E3 and termed QDQ, as a state intermediate between hadrons and 
2SC, but an identification with the 2SC may also be possible when the masses go 
to zero at high density. Although strictly speaking this QDQ and 2SC with zero 
net strangeness would be ustable against weak interaction decays in neutron stars, 
the number of strange quarks that is produced in the mixture can be neglected 
in a first approximation. The finite mass m s ~ 150 MeV suppress the ocupation 
of the corresponding Fermi sea, while the occupation of the bosonic QDQ or 2SC 
ground state is strongly favoured energetically. We must note that a similar mass 

1 /3 

parametrization ( m s — m s o + C/ rig ) has been employed to describe strange quark 
matter within the QMDD model. Also note the slightly different parametrization 
than the one employed in Ref. El motivated by an intuitive construction of the mass 
formula. 

In view of the above, we shall proceed to model the quark - diquark plasma as 
a gas of u and d quarks, and (ud) diquarks. Non-perturbative effects are accounted 
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for using the quark mass-density-dependent model of confinement to derive the 
properties of the equation of state. 

The expression for the thermodynamic potential density of particles of mass mi, 
chemical potential fii and degeneracy factor gi in the large volume limit of a free 
gas is 

Q. = g f°° l n (l ± e (^-(P 2 +™?) 1/2 )/^))p 2 dp (3) 

B in the QMDD model 
The pressure is actually 



o 



Because of the dependence of the mass with 
the standard relation P = —CI is no longer valid, 
given by 



For the other thermodynanhcal quantities we write (see Ref.[12]) 

E = n + Y J ^n l -T— (5) 



an 
dm 



(6) 



T,m k 



with E the energy density and ^i-the particle number density of the i-species. 

As it has been stressed before cm, for the fermionic component the term ub ^p*- ^ L 
forces the pressure to be zero at finite baryon number, playing the role of the bag 
constant B of the MIT Bag Model. Needless to say, this is a property related to the 
mass parametrization only, and has no relation with the statistics of the particles 
themselves (Bose or Fermi). However, it happens that in the case of bosons this 
extra term does not contribute at zero temperature. Thus, as in the standard case, 
for massive bosons at T = the energy density reduces to E = J^. nimi. 

As we mentioned above we shall model the quark - diquark plasma at zero 
temperature as a mixture of a Fermi relativistic gas of u and d quarks, and a Bose 
condensate of (ud) diquarks. We neglect the contribution of antiparticles and that 
of a leptonic component such as electrons, muons or neutrinos. This mixture is then 
described by B0 

P = ^ fninixlG(xi) - 2J m i n i f(x i ) (7) 

i—u.d i—u,d 

Diquarks do not contribute to the pressure in this approximation, since they are 
all in the ground state, but do contribute to the energy density 

E = X! m i n i F ( x i) + m D n D . (8) 

i— u : d 

The number density of quarks u and d is given by 
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where the functions F, G and / are given by 



BE 



G(xi) = I 



Xi(xj + 1) 1/2 (2a;- + 1) - argsinh(xj 



x l {x 2 i + l) 1/2 {2x 2 - 3) + 3 argsinh(xj) 



3 riB drrti 
2 m,i driB 



Xi(x 2 + l) 1 / 2 - argsinh(a;i) 



14 



14 



14 



(9) 



(10) 



(11) 



(12) 



being x = pp/rrii = (/j, 2 — mf) 1 / 2 /m,, \Xi the chemical potential and rrii the mass of 

the particle, which in our case depends on the baryon density. 

These equations must be complemented with the following conditions: 

1) Chemical equilibrium: D ^ u + d requires [id = A*« + fJ-d- If we use the 

fact that at T = the condensed diquarks all have zero kinetic energy we have 

(p D = m D ) 

m D =fiu + fJ-d- (13) 

Then, using the fact that, /i = (x 2 + l) 1 ' 2 m, m u = m^, and mry = hido + 2m M , 
we can write the chemical equilibrium as 



moo 



+ 2 = (xi + iy/ 2 + (x 2 + iy/ 2 . (u) 

l/3ra<z = 0, therefore, 



2) Electrical charge neutrality implies l/3nn + 2/3n 
the diquark number density is given by 



no = n d ~ 2n y 



;ml(x 3 d -2x 3 u ). 



3) Baryon number is given by ns — 2/3nr) + l/3(n u + rid), thus 
n B = n d - n u = \ml{x 3 d - x 3 u ). 



(15) 



(16) 



The equation of state can be solved numerically for each ub once the parameters 
(C, moo) are given (see Figs 1 and 2). 

As the baryon density increases the EOS tends to a pure "quark d-diquark" 
composition since the number density of quarks u goes to zero and the abundances 
of quarks d and diquarks D become equal. The same behaviour arises for smaller 
enough values of the mass mco (see below). In these cases the EOS takes a very 
simple form: x u — 0, n u — 0, ub — nu = n d — (m^a;^)/7r 2 . 

As expected from simple physical considerations, we find that the EOS of the 
quark-diquark is quite soft, especially at relatively high densities, because of diquark 
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condensation. The EOS becomes somewhat stiffer at low densities when compared 
to the MIT-inspired EOS developed in Ref.[18], at least for a large region of the 
parameter space. The same effect is found if we compare the MIT and the QMDD 
model EOS for strange quark matter uB. 

3. Stability of the quark - diquark matter 

In order to determine the allowed values for C and moo we impose the following 
stability criteria for quark-diquark matter. 

1) Instability of two flavor quark matter: The energy per baryon of two flavor 
quark matter (at P = and T = 0) must be higher than 939 MeV (the neutron 
mass) . Assuming that the current mass of quarks u and d is zero this results in the 
condition C > (159.3 MeV) 2 . 

The energy per baryon of pure udd matter at P = and T = is, for a given C 



E 
n B 



= m u [F{x u o) + 2F{x do )} (17) 

idd 



where m u is the mass of quark u in two flavor quark matter at zero P and T, 

c c cvv/ 3 



1/3 1/3 1/2 \^°> 



i-p nint 

300. 



and Xdo — 2 l ^ s x u o, where x u q = 1.13173096 is the value that produces a zer 
pressure, P u dd = 0. This criterion has already been used in previous papers E 

2) Relative stability of quark-diquark matter: The energy per baryon of quark - 
diquark matter (hereafter ud-D matter) must be lower than that of pure two flavor 
quark matter (at P = 0, T = 0). 

The zero pressure condition for quark - diquark matter at T — reads 

x 5 u G(x u ) - x 3 J(x u ) + x 5 d G(x d ) - x 3 J(x d ) = 0. (19) 
The energy per baryon is given by 

E x 3 u F(x u )+x d F(x d ) + (^ + 2)(x d -2xl) 
— = m u 3 _ 3 (20) 

U B ud-D X d ~ X u 

Note that the mass m u of quark u in ud — D matter is in fact a function of C, x u 
and x d given by m u — C/n 1 ^ 3 = C'/(n d — n.^) 1 / 3 , from which we derive 



(iS-i»)v«- 



Also, the factor (™_Do/ m u + 2) in Eq. J2Q| ) is a function of x u and x d given by Eq. 

(ED- 

The fundamental relation for this stability requirement is 
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E E 

n B U d—D n B 

therefore, from eqs. ( p0| ) and ( fl7j ) we obtain 

xlF{x u ) + x d F(x d ) + (^ + 2)(x d - 2x1) 



(x d - 4)V6 



(22) 

udd 



< n (23) 



where 



[F(x uQ ) + 2F(x d0 )} 
kq = — 2 . (24) 

x u0 

From eqs. ( p3| ) and (|l^) we find numerically the set of values x u \ and Xdi that 
satisfy the imposed stability conditions. 

Using chemical equilibrium alone (eq. (fbi])) we find a very simple relation be- 
tween ttidq and C 

m D0 = KC 1/2 (25) 

with 

g _ 7r i/ 3 (^i + l) 1/2 + (^i + l) 1/2 -2 f2gl 

C T 3 _ „3 U/6 ■ ^1 

This is an "isochemical" straight line along which the energy per baryon increases 
with C. The condition {E/n B )\ud-D < {E/n B )\udd is fullfilled for K < 1.723 (the 
upper limit is given by the full straight line in Figure 3). Also, we find that for 
K < 0.8417 there are no quarks u in the mixture and the composition is always a 
pure "quark d-diquark" mixture (see dashed line in Figure 3). 

3) Absolute stability condition: we have also looked for a region inside which 
quark - diquark matter could have even a lower energy per baryon than the neutron 
mass 



E 
n B 



< 939MeV, (27) 

ud-D 



being though absolutely stable. In the (C 1 / 2 ,mr)o) parameter space this region 
is " triangle- like" . The same shape has been obtained for strange quark matter in 
the QMDD model On the right side boundary of this " triangle- like" region the 
energy per baryon of quark-diquark matter at zero pressure is exactly 939 MeV. 
The region is shown in Fig. 3 and labeled with A and A . It is interesting to note 
that negative values of moo are allowed (see Fig. 3) and so the diquark may be 
bound in this model, that is, its mass will be smaller than the sum of the masses of 
its constituent particles. 



4. Quark-diquark stars 
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The exploration of the mass-radius relation for compact stars is in fact a potentially 
powerful tool in testing the existence of different phases of matter inside NSs, and 
their composition. On the other hand the observational data of the mass-radius 
relation imposes severe constaints on high density equations of state. 

Almost all of the meassured masses of " neutron stars" lie within a narrow range 
around 1.4 Mq. This mass scale may be due to the fact that neutron stars are formed 
in the gravitational collapse of supernovae and come just from the iron cores of core 
collapse supernovae. However, we may wonder whether smaller neutron stars could 
exist and whether there exist some mechanism capable of create them because these 
are not expected in core collapse models and hold the potential for discriminating 
among possible equations of state. This question becomes more significant in view 
of the recent claim of very low- mass/low-radii objects. At least two objects, Her 
X-l and RX J1856-37 are candidates to high compactness. They have been claimed 
to have radii ~ 7km and a masses around IMq 

If quark- diquark matter has a lower energy per baryon than normal nuclear 
matter we may expect the existence of stable stars made up entirely by this phase. 
By solving the Tolman-Oppcnheimcr-Volkoff equations of stellar structure, we ex- 
plored the mass- radius relation for quark-diquark stars assuming a parameter set 
inside the window of absolute stability. 

In Figs. 4 and 5 we see that for particular values of the parameters (C 1 / 2 , moo) 
the obtained values of M and R are compatible with present determinations of the 
above sources. We also find that the effect of increasing moo for a fixed C is to 
lower both the maximum radius and the maximum mass of the stellar configurations. 
Similarly, increasing C at fixed moo allows configurations with smaller maximum 
radius and mass. The (C 1 / 2 ,to_do) values selected for the sequences shown in Figs. 
4 and 5 match naturally the observed compactness of Her X-l and RX J1856-37. 

5. Conclusions 

We derived an equation of state for quark diquark-matter in a self-consistent phe- 
nomenological model where the mass of quarks and diquarks is parametrized with 
the baryon number density. The confining forces are included via a dependence 
of the mass with n B ^ 3 . The model depends on two free parameters (C 1 / 2 , moo) 
whose values can be limited invoking stability criteria and allows the existence of 
a stable state of the quark-diquark gas at zero pressure, for densities of the order 
of few times the nuclear saturation density (see Figs 1 and 2). The free parameter 
moo is introduced in order to describe in a simple form the unknown binding of 
the diquark. Negative values of moo should mimic the behaviour of bound color 
superconducting quark matter while positive values are considered in view of recent 
results i claiming a diquark mass of ~ 700 MeV which suggest that this particle 
is unbound. The mixture contains a Bose condensate of (ud) diquarks so that, the 
resulting equation of state is in fact very soft, as it can be seen in Figure 1 where 
we have depicted P versus E for different parameters . The here described quark- 
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diquark phase has always a non-zero fermionic pressure, since due to the excess 
of quarks d over quarks u there will always exist a free Fermi gas of quarks d to 
satisfy the equilibrium conditions. When compared with the equation of state for 
quark-diquark matter presented in E3 we found that independently of the value of 
the parameters, the latter is considerably stiffer at low densities. However, in some 
cases (c,d, e and f of Fig. 1) there is a limiting density beyond which the new equa- 
tions of state are softer than the previous one. This is relevant for compact star 
structure because it allows the existence of stable compact stars with even smaller 
radius than those calculated in previous models Ea. It has been argued that the 
few percent effects (O (A//i) 2 ) from the condensation are impossible to disentangle 
from other uncertainties, but the fact is that, for example, the whole strange matter 
theory is also subject to this caveat, and in this respect the models here presented 
open a qualitatively similar scenario for compact stars. 

We have found the regions in the parameter space (C 1 ' 2 ,m j Do) where quark- 
diquark matter in bulk is stable against decaying to pure two flavor quark matter 
and to a neutron gas (see Fig. 3). The region of absolute stability allows the 
existence of pure quark-diquark stars. It is important to note that the existence of 
self-bound compact stars is possible for both bound and unbound diquarks, since 
the assumed bosonic character of the diquark favors the transition of free quarks at 
the Fermi surface to the bosonic ground state and ensures a global energy decrease, 
more specifically a lower value of the energy per baryon, for the appropiate choice 
of parameters. That is, since typical Fermi energies in a two flavor quark mixture 
is ~ 300 MeV, the energy gained by condensing bosons and sending them to the 
ground state is in fact more relevant for stellar structure than the hypothctically 
gained energy of the pair gap. Moreover, even in the case of unbound diquarks 
a lower energy state is reached by the system as a whole. Therefore, although 
the parameter that mimics the pair gap may be not only negative but also zero 
or positive, energy is gained in all cases due to the assumed bosonic character of 
the superconducting pair. A slightly different parametrization used in @ produced 
different results. In particular the resulting stellar models were not so compact, a 
feature due to a larger stifness of that equation of state for all the parameter range. 
This gauges the sensitivity of the models to the assumptions on the basic physics. 

Finally, we must note that, strictly speaking, quark-diquark matter is unstable 
against weak interactions producing strangeness. In this sense, the present study 
may be considered as a simplification analogous to pure neutron matter models of 
neutron stars. A more general study should include the possibility of strangeness 
production, the appeareance of leptons in the mixture, as well as 2SC and color 
flavor locking (CFL). We can expect that the opening of the weak channel will 
produce an even softer equation of state and so higher compactness of these stars. 
These improvements will be presented in future work. 
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Fig. 1. Equation of state of quark diquark matter in the quark mass density dependent model 
for different set of parameters (C 1 ^ 2 , mjjo)- The curves are labeled with: a = (160,-100), b 
= (160,0), c = (200,-100), d = (160,210), e = (185,100) and f = (185,130). For the sake of 
comparison we also include the ultrarelativistic EOS P = E/3 and the EOS presented in Rcf. [18] 
(labeled as HP). 



Fig. 2. The particle number densities of u and d quarks and diquarks for the parameters (160, 200) 
(label A) and (185, 130) (label B) . At high baryon number densities the number of quarks u goes 
to zero and the number of quarks d and diquarks are equal. We also include the case of pure 
quark d - diquark matter tib — = nrj corresponding to any set of parameters in the region A' 
of Figure 3. 



Fig. 3. Different stability regions in the parameter space of the model. Inside regions A and A 
quark-diquark matter has a lower energy per baryon than a neutron being so absolutely stable. 
For parameters inside the region A there are no u quarks in the mixture, and matter is always 
composed by quarks d and diquarks. Inside the region B below the full straight line the mixture 
has a lower energy than pure two flavor quark matter composed only by quarks u and d. 



Fig. 4. The mass radius relation for quark-diquark stars. The curve labeled with HP corresponds 
to the EOS developed in E3- The other curves correspond to the following set of the parameters 
{C 1 ^ 2 ,m D0 ) : a= (185,130), b=(185, 100), c=(180, 100), d=(160, 210). Models overlap the range 
of the M-R measurements of Her X-l and RX J1856-37. 



Fig. 5. The same as the previous figure but for A= (210,-100), B= (200,-100), C= (160,0), 
D=(160,-100). 
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